Dephasing by a nonstationary classical intermittent noise 
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We consider a new phenomenological model for a classical intermittent noise and study 

its effects on the dephasing of a two-level system. Within this model, the evolution of the relative 
phase between the |±) states is described as a continuous time random walk (CTRW). Using renewal 
theory, we find exact expressions for the dephasing factor and identify the physically relevant various 
regimes in terms of the coupling to the noise. In particular, we point out the consequences of the 
non-stationarity and pronounced non-Gaussian features of this noise, including some new anomalous 
and aging dephasing scenarii. 
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I. INTRODUCTION 

Recent experimental progress in the study of solid- 
state quantum bits (Josephson qubits)i has stressed the 
importance of low-frequency noise in the dephasing or 
decoherence of these two-level systemsS*^. It now ap- 
pears that the coupling to low-frequency noise is the main 
limitation in obtaining long lived phase coherent states 
of qubits necessary for quantum computation. However, 
a complete understanding of the microscopic origin of 
1/ / noise in solid state physics is not available yel£ and 
therefore, theoretical studies of the dephasing by such a 
noise are based on phenomenological models. In the spin- 
boson model, the environment of the qubit is modeled by 
a set of harmonic oscillators, with an adequate frequency 
spectrum-Si Another commonly used model for a low- 
frequency noise consists in considering the contributions 
from many independent bistable fluctuators^. In the 
semi-classical limit, the noise from each fluctuator is ap- 
proximated by a telegraph noise of characteristic switch- 
ing rate 7. For a broad distribution ~ I/7 of switching 
rates 7, a 1// spectrum is recovered when summing con- 
tributions of all fluctuators. Such a model is based on 
observations of telegraph like fluctuations in nano-scale 
devicesifliii, but a precise characterization and justifica- 
tion of the broad distribution of switching rates is still 
lacking although the localization of these fluctuatorsiSii^ 
as well as their collective or individual nature^ have been 
investigated for a long time. 

In this paper, we consider a new phenomenological 
model for the classical low-frequency noise. This model 
can be viewed as the intermittent limit of the sum of 
telegraphic signals. In this limit, the duration of each 
plateau of the telegraphic signal is assumed to be much 
shorter than the waiting time between plateaus^. A 1// 
power spectrum for the intermittent noise is then recov- 
ered for a distribution of waiting times r behaving as 
r~ 2 for large times. Because the average waiting time is 
infinite, no time scale characterizes the evolution of the 
noise which is nonstationary. The purpose of this article 
is to study the effects of such a low-frequency intermit- 
tent noise on the dephasing of a two-level system in order 
to identify possible signatures of intermittence. 



As we will show, in this model the relative phase <& 
between the states of the qubit performs a continuous 
time random wallsi^*^ (CTRW) as time goes on. Such 
a CTRW was considered in the context of 1// current 
noise by Tunalejiii, extending the previous work of Mon- 
troll and Schcr on electronic transport^. However, in 
the present paper it is the integral of noise, and not 
the noise itself, which performs a CTRW. Moreover to 
our knowledge, the precise consequences of CTRWs non- 
stationarity on dephasing have not been studied. On 
the other hand aging CTRW were previously considered 
in the context of trap models in glassy materials^ and 
in the study of fluorescence of single nanocrystals22i2I. 
Technically, the dephasing factor that we will consider 
corresponds to the average Fourier transform of the po- 
sitional correlation function of the random walk. Some 
of the asymptotic behaviors of this correlation function 
were already obtained in ref^J However, in the present 
paper we will extend these results to all possible regimes 
and we will present all of these results in a unified frame- 
work. The use of renewal theorjiSS greatly enlightens 
the origin of non stationarity and enables us to interpret 
some features of the dephasing scenarii. 



This paper is organized as follows : in section [H] we 
present our model for the noise and define the quantity 
of interest, i.e. the dephasing factor of a two-level system 
coupled to this noise. In section IIIII the exact expres- 
sion for the single Laplace transform of the dephasing 
factor will be derived and, from this result, the physi- 
cally relevant weak and strong coupling regimes are iden- 
tified. Moreover, we clarify the origin of non-stationarity 
and show the relation of our problem to renewal theory. 
For completeness and pedagogy, the effects of standard 
anomalous diffusion of the phase and of randomness of 
waiting times on dephasing are compared showing the 
importance of intermittence in the non-stationarity prop- 
erties of the dephasing scenarii. In sections HVl and Ivl we 
present a complete study of the behavior of the dephasing 
factor respectively for a noise with a vanishing average 
amplitude (symmetric noise) and with a finite average 
one (asymmetric noise). The general discussion of the 
results is postponed to section IVT1 
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II. THE MODEL 
A. Pure dephasing by an intermittent noise 

In this paper, we consider a quantum bit defined as a 
two-level system with controllable energy difference Tuoq 
and tunneling amplitude A between the two states |— ) 
and |+) (eigenstates of a z ). The effect of the environ- 
ment on this two-level system will be accounted for by a 
fluctuating shift TiX of the energy difference Tiuj$. Thus 
the Hamiltonian describing this model is written as 

H=^(uj cT z + A<j x -Xa z ) , (1) 

In this paper, we will mainly focus on the case of pure 
dephasing (A = 0). However, as explained below in sec- 
tion [^21 our discussion will also apply to other operat- 
ing points (A 0), including the special points where 
a careful choice of control parameters considerably lower 
the qubit sensitivity to low frequency noised. 

Here, we will focus on the effects of a low-frequency 
classical noise on the qubit. The noise is represented 
by a classical stochastic function corresponding to the 
fluctuations of the noise in a given sample. Within this 
statistical approach, we focus on the statistical properties 
(e.g. the average) of physical quantities associated with 
the qubit such as the so-called (average) dephasing factor. 
As we shall see now, its meaning can be understood by 
considering a typical Ramsey (interference) experiment 
on the qubi1p2i. 

In such an experiment, the qubit is prepared at initial 
time t p in a superposition of the eigenstates of a z , e.g 
l+) = (IT) + U))/\/2- N ote that throughout this paper, 
t = will correspond to the origin of time for the noise 
(e.g the time at which the sample reached the experi- 
ment's temperature). At some later time t p + r exp > t p , 
we consider the projection of the evolved qubit state on 
In the mean time, the state has evolved under Hamil- 
tonian (A = 0) and both states ||) and ||) have accu- 
mulated a random relative phase &(t pi T exp ) defined by 

$(t pi T exp ) = / X(t)dt. (2) 

Jt p 

For a given accumulated phase $ = &(t p , T exp ), the quan- 
tum probability P$^ T(!xp (||)) to find the qubit in state |+) 
at time t p + r exp is given by 

- P *:T exp (|+)) = - [1 + COS (uJoTexp - $(t p ,T exp ))} (3) 

Note that in a given sample, ■P$,r e ™(|+)) oscillates be- 
tween and 1 as a function of r exp (although possibly 
nonperiodically). The experimental determination of the 
probability for finding the qubit in the |+) state at time 
t p + t usually requires many experimental runs of same 
duration r exp . The phase fluctuations between different 



runs induce an attenuation of the amplitudes of these os- 
cillations (analogously to destructive interference effects 
in optics). Using Bayes theorem, the corresponding sta- 
tistical frequency to find the qubit in state |+) after a 
duration r exp is given by the probability 

Pt p ,r exp (\+)) = J d ® *V«p(l+»P(* = ®(t p ,T exp )) 

= \{l + $l{D tp (T exp )e- i »^)) (4) 

In this expression, the decay rate of these oscillations is 
encoded in the dephasing factor D tp (T exp ) defined as 

Dt p (T eX p) = exp(i$(i p , T ex p))- (5) 

In this formula (and only here), the overline denotes an 
average over all possible configurations of noise X(t) dur- 
ing the experiment. 

Note that in deriving eq.Q), statistical independence 
of the phases $ between different runs has been assumed. 
This is not necessarily true for successive runs in a given 
sample as correlations of the noise might lead to a depen- 
dence of the distribution of the phase &(t p , r exp ) on the 
starting date t p of the run. Hence throughout this pa- 
per, for self-consistency, we will keep track of this effect 
through a possible t p dependence of the dephasing fac- 
tor D tp (T exp ). Its possible implications will be discussed 
together with our results in section IvTl 

B. A model for classical intermittent noise 

In several experimental situations, the low-frequency 
noise acting on the qubit is supposed to be due to con- 
tributions from background charges in the substrate^. 
When the dephasing is dominated by the low-frequency 
fluctuators, a semi-classical approach, in which the noise 
is modeled by a classical field, appears sufficien1i£*2*2^. In 
this case, the noise is described by a Dutta-Horn modelSi. 
In its simplest form, the potential X(t) is written as 
the sum of the contributions of many telegraphic sig- 
nals, each with a characteristic switching rate 7 between 
the up and down states (see figure^). For switching 
rates distributed according to an algebraic distribution 
75(7) ~ I/7, the power spectrum of the corresponding 
noise has a 1// low-frequency behavior. 

In this paper, motivated by several noise 
signaturesiii^, we will focus on the intermittent 
limit for such a Dutta-Horn model, and propose a 
phenomcnological spike field to study its effects on qubit 
dephasing. By intermittence, we mean for the case of 
the telegraphic noise that the switching rate 7+ from the 
up to the down states is much larger than the switching 
rate 7_ from the down to the up states (or vice- versa). 
In this limit, the total noise reduces to a collection of 
well defined events, separated by waiting times Tj (sec 
figure 2J)). If considered on times much longer than the 
typical duration tq of these events (or at frequencies 
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iM 1] - 1/7 W 



(2) 



(2) 



1/7- > 1/74 




(a) (b) 

FIG. 1: Representation of a low- frequency noise as a sum of contributions from telegraphic signals (a). In this first case, the 
switching rates for the up (7+) and down (7-) states are comparable, and a I// spectrum is recovered for a distribution of 
switching rates ~ 1/7. The intermittent limit (b) corresponds to the limit where the noise stays in the down states most of the 
time (7+ 3> 7-). In this paper, we will approximate this intermittent noise by a spike field. For this intermittent noise, a 1/f 
spectrum implies a non stationarity whose consequences on dephasing are studied in this paper. 



smaller than 1/to), we can approximate this noise by a 
spike field consisting of a succession of delta functions of 
weight Xi corresponding to the integral over time of the 
corresponding events of the intermittent field (FigGJj). 

More precisely, denoting by t — the origin of time, 
and by tj the successive waiting times between the spikes 
(or events), we know that the nth spike occurs at time 
t n = y\-| Tj. The value of the stochastic intermittent 
classical noise X(t) (see Fig. 0) is then 



(6) 



x(t) 



FIG. 2: Representation of the random spike field used to 
model the intermittent low-frequency noise in this work. This 
noise is described by the distributions of the phase pulses Xi 
and of the waiting time intervals r% between the spikes. 



In the following, we will consider the dephasing pro- 
duced by this spike field. We expect any short time de- 
tails like the specific shape of the real pulses to be ir- 
relevant in the limit of typical dephasing time long com- 
pared to tq. Within this approximation, a noise signal is 



fully characterized by the collection of waiting times 73 
and pulse amplitudes Xi that occur as time goes on. We 
will assume these two quantities to be independent from 
each other, and completely uncorrelated in time. We will 
then characterize such a noise solely by two independent 
probability distributions ipi T ) an< l p( x ) f° r the Tj and Xi 
respectively. 

The distribution of pulse amplitudes p(x) will be as- 
sumed to have at least its two first moments finite, and 
denoted in the following by 



xp(x) dx , 



g = V x 



x 2 p(x) dx 



1/2 



(7a) 
(7b) 



We will consider separately the case of zero average (sym- 
metric noise) and of non-zero average (asymmetric noise) 
since the latter induces specific features of the dephasing 
factor, discussed in section IV! 

We will consider algebraic distribution of waiting times 
ip(r). parametrized by a single parameter ji : 



TQ 



TO 



TQ 



1+M 



T > 0. 



(8) 



As we will show in section Hi Dl a 1// M power spectrum 
for this intermittent spike field follows naturally from 
such a choice for VK 7 ")- As wc shall see in this paper, 
the above algebraic distribution of waiting times allows 
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to correctly capture the essential features of the dcphas- 
ing scenarii generated by an intermittent noise. The first 
and second moments of ip are finite respectively for // > 1 
and [i > 2 and are given by 

to = 7^7 ; < r2 > = 7 ^—^A W 

fi-1 (ix-l)(n-2) 

Let us note that the above defined pulse noise contains 
two independent potential sources of dephasing: The ran- 
domness of the pulse weights and the randomness the 
waiting times. Their respective effects will be compared 
in section ITTT1 

C. Decoherence at optimal points 

Before turning to the detailed study of pure dephas- 
ing (A — in QJ), let us mention that our discussion 
can be easily extended to the study of dephasing in the 
presence of a transverse coupling in in particular at 
the so-called optimal points. They correspond to con- 
figurations where the fluctuations of the effective qubit 
level splitting are only quadratic in the noise amplitude. 
The qubit can be operated at these optimal points by 
a careful choice of the control parameters ujq and A of 
the qubit and then, the influence of low frequency noise 
can be reduced considerably^. For the Hamiltonian con- 
sidered in the present work Q such an optimal point 
is reached for transverse coupling to the noise (cjq = 
and A / 0). In this case and assuming the amplitude 
of the noise to be small compared to the control pa- 
ra mctcr A, th e effective qubit level splitting is given by 
V / A 2 +X{t) 2 fx A + X(t) 2 /(2A). Hence, the dephasing 
effect of a linear transverse noise can be accounted for 
using an effective quadratic longitudinal noise. The cor- 
responding dephasing factor is then given by (J5J and 
with the replacement X — > X 2 /(2A). 

In addition, the transverse noise at an optimal point 
induces transitions between the eigenstates of the qubit, 
i.e. it leads to relaxation. The Fermi Golden Rule relax- 
ation rate IV, which is used for estimating the relaxation 
contribution to dephasing, involves the power spectrum 
of the noise at the resonance frequency of the qubiti The 
total dephasing rate is obtained by summing the contri- 
bution of relaxation given by T r /2 and the contribution 
of pure dephasing due to the above effective longitudinal 
noise. 

In general the statistics of X and X 2 are very different 
and a special treatment is needed to derive the dephasing 
factor at optimal points^. However, for the noise con- 
sidered in the present work, the effective quadratically 
coupled longitudinal noise can be viewed as en effective 
linearly coupled noise of the same type but with renor- 
malized parameters. The rcnormalizcd distribution of 
the spike intensities is now given by (x > 0): 




where To is a microscopic time scale needed to regularize 
the square of delta functions. Therefore, our results for 
the longitudinal noise presented in section III Al can also 
be used to describe the effect of a transverse noise. 



D. Spectral properties of the intermittent noise 

1. One and two point correlation functions 

To make contact with other descriptions of low- 
frequency noise, we will determine the behavior of the 
two time correlation function of our noise, or cquiva- 
lently of its spectral density. However, as we will see, 
this spectral density is far from enough to characterize 
the statistics of the noise for small /i, in particular due to 
its non stationarity. We will consider for simplicity the 
case of a non zero average h = x. 

a. Time dependent average Let us consider the av- 
erage of the noise amplitude X(t). In our case, it reduces 
to two independent averages: over the amplitudes Xi of 
the spikes and over the waiting times t.; between them. 
Noting that X(t) vanishes except if a spike occurs at time 
t, we can express its average in terms of the average den- 
sity S(t) of pulses at time t also called the sprinkling time 
distribution in ref. 03 : 



(X(t))=hS(t). (11) 

Using the expressions <|A4|1 , IjAlOD for S (t) (see appendix 
|A"|) . we obtain for the average of X : 

Vm = ^r = -{fM-l) for M >2, (12a) 

V / T 

(12b) 

jm- S -^~(^ for „<1. (12c) 

7T To V t / 

Hence the non-stationarity of the noise manifests itself 
already in the time dependence of this single time aver- 
age. While it is only a subleading algebraic correction 
for 1 < n < 2, this time dependence becomes dominant 
for < fi < 1. 

b. Two time function Following the same lines of 
reasoning, we can derive the expression of the two time 
correlation functions (with t > 0) : 

(X(t p ) X(t p + t)) = h 2 S(t p )S(t) . (13) 

The first factor S(t p ) corresponds to the probability that 
a spike occurs at time t pi while the second factor S(t) 
reads the probability of having a pulse at time t p + t, 
knowing that there was one at t p . This reflects the reini- 
tialization of the noise once a spike has occurred at time 
t p . From (|11I13|I . we obtain the connected two points 



functions: 



C(t p ,t) = (X(t p )X{t p + t)) c 

= h 2 S(t p )(S(t)-S(t p + t)) 

2. 1/f noise spectrum 



(14) 



We will define the spectral density of the noise as the 
Fourier transform of the connected two points functions 
restricted to t > : 



Sx(t P , 



C{t p ,t) cos (wt)dt. (15) 



The correlation function C{t p ,t) generically depends on 
both times t p and t p + t thus showing that in general X 
is not a stationary process. However, to extract the low- 
frequency behavior of Sx(t p ,u>), it appears sufficient to 
consider the quasi stationary regime \t\ <C t p which cor- 
responds to experiment durations much smaller than the 
age of the noise. In this regime, the connected correlation 
function l|14|) reduces for t > to 



C(t p ,t)~h 2 S(t p ) (S(t)-S(t p )) 



(16) 



The associated effective power spectrum is defined for 
frequencies lo large compared to l/t p and reads: 

S x {t P ,oj) ~ 2h 2 S(t p )$l(L[S}(-iu})) (17) 

where L[S] denotes the Laplace transform of S(t). Note 
that in this quasi-stationary regime, the nonstationarity 
of the noise manifests itself only through the t p dependent 
amplitude S(t p ). Using explicit expressions for L[S] (see 
appendix 0), we obtain the effective power spectra: 



Sx(t p ,u)~ h 2 S(t p ) 
for < fx < 1, and 



cos (7T/X/2) 



r(i - m) 



(wT )~ 



(18) 



S x (t p ,(jj) ^ 

h 2 s(t p ) si n(^fVy»-ir-i 



("<r» 



fi — 2 



(19) 



for the intermediate class 1 < it < 2. The common 1/lo 
dependence of the spectral density is recovered in the 
limit fx — > 1. More precisely, the Laplace transform of S 
for /j, = 1, obtained in appendix lAl gives a logarithmic 
correction to a 1/u effective power spectrum: 



S x {t p ,uj) 



nh 2 S(t p ) 
lotq (log (wt )) 2 



(20) 



Let us stress finally that this effective power spectrum, 
although useful to compare our approach with other noise 
models, is not sufficient to characterize the statistical 
properties of the spike field noise relevant for dephas- 
ing. As we shall see in this paper, it does not account 
precisely for the non-stationarity of the dephasing factor. 
Besides this, non Gaussian properties of the noise have 
strong effects on the dephasing factor in many regimes. 
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FIG. 3: A configuration of the noise X(t p + r) in our model (as 
a function of r) between t p = 10 5 to and t p + r = 1.1 10 5 ro. In 
this figure, h = and the waiting times are distributed with 
an algebraic distribution P(ri) ~ t~ 21 . The bottom part of 
the figure shows the corresponding continuous time random 
walk of the accumulated phase of the TLS between t p and 
tp + T. 



III. DEPHASING, CONTINUOUS TIME 
RANDOM WALK OF THE PHASE AND 
RENEWAL THEORY 

A. Continuous time random walk of the phase 

Having defined our model for the intermittent noise 
field X(t), we will now study its effects on the dephasing 
of the two-level system, characterized by the dephasing 
factor ©. Note that a priori this dephasing has several 
origins : the randomness of the pulses amplitudes Xi , and 
the randomness of the waiting times Tj. These two de- 
phasing sources are assumed to be independent from each 
other in this work. For clarity and pedagogical reasons, 
we first start by considering the effect of randomness of 
the pulse amplitudes before turning to the the effect of 
waiting time randomness. 



I. Dephasing by random phase pulses 
the phase 



the random walk of 



Here, we assume that all waiting times are equal 
Tj = t- t — = f. In this case, the phase (fr(t pi t) ac- 
cumulated between t p and t p + 1 performs a usual ran- 
dom walk characterized by the distribution p(x) of phase 
pulses : at the dates t n = nf, <f>(t p ,t) is increased by a 
random value Xi. In the limit \x\ -C 1, the phase slowly 
diffuses and dephasing is achieved only after a large num- 
ber of pulses n = t n /f. Then, the distribution of the 
phase <I>(tp, t n ) can be well approximated by a Gaussian 
distribution (apart from some irrelevant tails) and the 
dephasing factor is easily computed. It decays exponcn- 
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tially with characteristic time = g i_ h 2 ■ 

Note that this diffusive regime can also be studied 
when the distribution p{x) lacks a finite second or even 
first moment. An anomalous diffusion of the phase is 
expected^ 9 .. Let us assume that a;, = /; + n£i where 
the probability distribution of £, has zero average and 
belongs to the attraction bassin of the stable law L v ^ 
characterized by the exponent < v < 2 of the algebraic 
tails for large values of £j and the asymmetry parame- 
ter |/3| < 1. Then, according to the generalized central 
limit theorem, the accumulated phase after n pulses is 
= nh + k n 1 / v £ where £ is distributed according to 
ii/,/3- Consequently, in the diffusive limit (n <C 1), the 
dephasing factor is the product of the homogenous phase 
e »ht/r £jj e characteristic function of L v ,p evaluated for 
k ~ git/f) 1 ^. For v < 2 and v ^ 1, this leads to: 

£) t (Jj = e ~CK, v t/T e i(h-Cl3ta,n(-Kv/2)K, v )t/f ^ ^l) 

where C is a numerical constant which can be absorbed 
in a rescaling of k. Thus, in this case, the decay is still 
exponential and stationary although non Gaussian fea- 
tures of the noise lead to an anomalous dependance of 
the dephasing time on the coupling constant k that char- 
acterizes the scale of the fluctuations of phase pulses^. 

2. Dephasing by random waiting times and continuous time 
random walk 

Let us now turn to the situation where the phase pulses 
happen at random times. In this case, the accumulated 
relative phase $(t p ,r) (see cq.(j2J|) does not perform a 
random walk as r increases, but rather a continuous time 
random walki^i^ (CTRW) on the unit circle. In other 
words, after some random waiting time Tj, ^>(t p ,r) is in- 
cremented by a random value Xi (see figure [3J. Thus, on 
a technical level, the dephasing properties of the two level 
system are now related to some correlation function of the 
corresponding CTRW. The corresponding dephasing fac- 
tor can differ from results obtained in the previous section 
due to the additional source of dephasing given by the 
randomness of waiting times. To illustrate this point, let 
us consider the case where all phase pulses have the same 
intensity p(x) = 6(x—h). The accumulated phase after N 
events is exactly Nh. Thus, dephasing comes only from 
the randomness of the number N[t p , t p + t] of events oc- 
curring between t p and t p + t and the phase diffusion is 
governed by the probability distribution for N[t p , t p + t}. 
The time needed to obtain N events is the sum over the 
iV first waiting times after t p : = Ylj=i T j- At t p = 0, 
all TiS are distributed according to the same probability 
distribution and therefore, in the limit of large N, the 
generalized central limit theorem^ can be used (the case 
t p 7^ will be discussed below) . It provides the limit law 
for tM at large N which in turn determines the proba- 
bility law for the number of events. According to this 
theorem, three classes must be considered depending on 
whether the moments @ are finite: 



1. The case /i > 2 where both (r) and (r 2 ) arc fi- 
nite. In this case, the probability distribution for 
the number of events is gaussian with a vanishing 
relative uncertainty. 

2. The case 1 < [i < 2 where the average of r is finite 
but the second moment diverges. 

3. The case \x < 1 where all moments diverge. 

The usual model for telegraphic noise assumes a Poisso- 
nian distribution for the number of events in a given time 
interval and corresponds to ip( T ) = je~ 1T . We will refer 
it as the Poissonian case and it belongs to the fi > 2 
class. 

Note that non trivial behavior is expected in the last 
two cases where ip has infinite first or second moments. In 
particular, for fi < 1, as ip(r) does not have any average, 
no time scale characterizes the evolution of the noise and, 
as we will see, nonstationarity follows. This case deserves 
a special attention as we showed that a 1/f spectrum is 
found precisely for fi —> 1. Before turning to the general 
formal expressions for the dephasing factor, we will focus 
on the origin of this nonstationarity. 

3. Origin of the nonstationarity in CTRW : the first 
waiting time distribution. 

Within our model, the waiting times between succes- 
sive pulses are chosen independently according to the 
distribution ip. Consequently, all the t p dependence of 
&(t p ,t) will come from the choice of t\ defined as the 
waiting time between t p and the first spike that follows 
t p (see figure 0J. Indeed, at time t p + n, the CTRW 



x(t) 







n 












— > 








to 




► 

t 



tp tp ~\~ Texp 



FIG. 4: Intermittent noise between t p and t p + t 

starts anew: T2 is chosen without any correlation to the 
history of the CTRW. Hence given the probability dis- 
tribution of Ti, we can forget about the history of the 
CTRW of the phase and describe its behavior starting 
at t p . This remark is at the core of the use of renewal 
theory. In the following, the probability distribution of 
Ti will be denoted by tpt„ and a priori, it may depend on 
t p . In fact, as we shall see later, its behavior can be quite 
counter- intuitive. 

First of all, note that a given t\ can be obtained from 
many different noise configurations that differ from the 
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time of the last event occurring before t p . Separating 
noise configurations (starting at t = 0) that have their 
first spike at time t p + n from the others leads to an 
integral equation that determines ip tp in terms of ip (S is 
determined from tp through an integral equation (|A1|) ~): 



<Mn) = i>{t P + n) + / dr ^(ri + T)S{t p - t) . (22) 
Jo 

The integral in the r.h.s. comes from noise configura- 
tions that have a spike between and t p . Equation 122|) 
is the starting point for deriving analytic results about 
ip t in appendix [5] an d [U] using Laplace transform tech- 
niques. Before computing exactly the dephasing factor, 
let us show some of the counter-intuitive properties of 
-0 tp (ri). In particular from (|22|l . we can derive the fol- 
lowing expression for the average of t\ valid for fi > 2: 



2(r 



2 r 



(23) 



This first term corresponds to the case of regularly spaced 
pulses averaged over the origin of times. The second term 
is the contribution of fluctuations. This result means that 
irregularities in the event spacings increase the average 
waiting time of the first event following t p . The r.h.s. of 
(I23|l does not depend on t p , as expected from the sta- 
tionarity of the CTRW for /i > 2 after a short transient 
regime at small t p . On the other hand, for /i < 2, eq. 
(|23[) is expected to break down since (r 2 ) diverges. This 
divergence signals that in some noise configurations Ti 
can become of the order of t p and, as a consequence, 
the average properties of the CTRW after t p , depend on 
this age of the noise. Indeed, we can show from i|22|) that 
(ri)^, t scales with the age as t^~^. Note that in the diffu- 
sion regime, many phase pulses are necessary to dephase 
the qubit. Therefore, we expect that, in this regime, the 
aging effect on ipt p (Ti) only brings weak corrections to 
the dephasing scenario as will be confirmed below by ex- 
act computations. However, we shall see in the following 
that the t p dependence of ip t for fi < 1 has much more 
spectacular consequences on the dephasing factor than in 
the 1 < [i < 2. 

We will now show that, knowing the Laplace transform 
of ipt v , an explicit expression for the Laplace transform 
of the average dephasing factor can be obtained. 



B. Exact dephasing via renewal theory 



All other histories have at least one event between t p and 



t p + t. Let us assume that it happens at t p 



t where 

t lies between and t. Then after this event, the noise 
starts anew. The jump itself contributes by e lx to the 
dephasing factor and the rest of the noise configuration 
contributes by D (t — r) (i.e. t p = in eq.JSJ). The 
probability that the first event after t p happens at time 
t p + t is nothing but fpt p (T). Hence the contributions 
to D tp (t) from all possible noise configurations take the 
form of the following renewal equation: 



D tp (t)=U Q (tp,t)+e" / ip tp ( T )D (t-T)dT. (25) 
Jo 

Note that the Laplace transform of this expression is very 
simple: 



L[D tp ] = Lpo] + (!-/) (XfofrtJ. L[D ]) 



(26) 



where / = 1 — e ix . Specializing t p = 0, one gets an ex- 
pression for D (t) that contains n (0, t). Since n (0, t) = 
J.°° ip(r) dr an explicit expression for the Laplace trans- 
form of Dq can be found: 



L[D ] 



1 1 - L&] 



s 1-(1 -/)!#]■ 
Plugging expression i|27|) into equation (|26fl gives: 



L[D tp ] 



1 



1 



1-(W)£[ 



(27) 



(28) 



This exact expression will be extensively used to derive 
both analytic expressions and numerical plots by Laplace 
inversion. Before proceeding along, let us express L[D t ] 
in terms of the density of events. It is then convenient 
to define St {t) = S(t p + t) (see appendix^}. From the 
expressions in this appendix of i[<S't p ] and L[S] in terms 
of L[tp] and L[ip tp ], we obtain 



L[Dt p ] = 



1 



1 - 



fL[S tp ] 
l + fL[S] 



(29) 



Equation l|28|) and 129fl enable analytical estimates of the 
dephasing factor in various limiting regimes. Equation 
(|29|l is useful in the limit of small coupling (i.e. /—►()) 
whereas eq. I|28f) will be useful in the opposite limit of a 
wide distribution of spikes' heights (i.e. /—»!). 



1. Dephasing factor 



2. Weak and strong coupling regimes 



Among all noise configurations that are to be taken 
into account, some of them (possibly very few) do not 
have any event between t p and t p + 1. Their total weight 
is given by Ho(t p , t) which is 



/+oo 



^t p (r) dr. 



(24) 



In the limit of very strong coupling (/ = 1), the 
phase spreading of a single spike is sufficient to dephase 
the qubit. In this case, ea. H25|) immediately leads to 
D tp (t) = Ho(t p , t). This means that the whole average is 
dominated by the rare noise configurations that do not 
evolve during the experiment. As we shall see later, this 
leads to a strong t p dependence of D tp (t) for ji < 2. 
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At much lower couplings (|/| <C 1), it is not a priori 
clear whether the renewal equation for D tp (t) induces a 
strong dependence of D tp (t) on t p . Pushing the above 
analysis forward amounts to compare the time depen- 
dence of the two terms in the r.h.s. of the renewal equa- 
tion l|25|) . As we will see now, it will provide a better 
understanding of the physics underlying the dephasing 
scenario. 

Let us first assume that U (t p ,t) decays much faster 
than D (t). Since ipt p ( T ) = —d T U (t p ,T), it also decays 
much faster than D . Consequently, we can approximate 
i>t P (r)D (t-T) by i) tp (T)D {t) in Eq.©. Then, after a 
short initial decay due to both D (t) and n (< p , t), D tp (t) 
decays as Do(t). Consequently, the dephasing time does 
not depend on t p . For instance, we expect this situation 
to occur at weak coupling when the average waiting time 
(r) is finite since the probability that no event occurs be- 
tween t p and t p + t vanishes quite fast when t increases. 
Note that in the limit of vanishing coupling, the spread- 
ing of the phase can become much slower than the de- 
cay of n (t p , t) and thus our starting point hypothesis is 
valid. 

The opposite case where the decay of Do(t) is much 
faster than the decay of Uo(t p ,t) can be discussed more 
conveniently by integrating the renewal equation (|25l) by 
parts. Considering again g 1 we get: 



D tp (t) ~ D {t) - / U (t p ,t-r)D' {T)dT. 



(30) 



Approximating Tlo(t p ,t — r) « IIo(tp,t) in this equation 
yields D tp (t ) = D a (t) + n (t p , t ) ( 1 - D a (t ) ) . As a conse- 
quence, once D (t) has decayed, the dephasing factor is 
almost equal to Ho{t p , t). This is the same behavior than 
in the very strong coupling regime / = 1, although here, 
we assumed |/| <C 1. This regime is expected to occur 
when the decaying time scales of Do(t) and of Ho(t p ,t) 
are comparable. This is obviously the case at very strong 
coupling but, surprisingly, as we shall see now it can also 
be obtained for |/| <C l! 

First of all, the above discussion shows that such a 
regime can only happen if the average waiting time is 
infinite, i.e. for fi < 1. In this case, Ho(t p ,t) can be 
evaluated analytically (see eq. (|D 1|> ^ : it is shown to be 
independent of To and to exhibit aging behavior (i.e. to 
depend only on t/t p ). Therefore, the decaying time scale 
of IIo(ip, t) scales as t p . Comparing this time scale with 
the dephasing time scale for |/| <C 1, leads to a t p de- 
pendent cross-over coupling constant f c {t p ). In the case 
of aging IIo, the crossover coupling decays to lower and 
lower values by increasing t p . 

For |/| > fc(t p ), the dephasing factor behaves as 
Ho(tp,t) and the dephasing time saturates func- 
tion of /. Such a saturation of the dephasing time as 
a function of the amplitude of the noise has already been 
discussed for a Poissonian fluctuates*^. In this case 
n (i p ,t) = n (0, t) = e~*/^ decays very fast, on a time 
scale (r). The cross-over between weak and strong cou- 
pling regime happens precisely at the point where the 



dephasing time assuming weak coupling (r) / g 2 is of the 
same order as this decay time, i.e. for g ~ 1. Note that in 
this case, as expected from our discussion, the crossover 
scale is independent of the age of the noise t p . We expect 
this reasoning to break down in our model because of the 
broad distribution of waiting times (jx < 2). Understand- 
ing the various dephasing scenarii and computing the t p 
dependence of the crossover coupling in our case requires 
the computation of iIo(t p , t) and of Do(t). These quanti- 
tative results will be presented in forthcoming sections. 

To summarize the above discussion, we have argued 
that the dephasing time is bounded by the typical decay 
time of both Ho(t p , t) and D (t). This suggests to distin- 
guish between two regimes: on one hand, a weak coupling 
regime for which Do(t) decays much slower than Tlo(t p , t) 
and for which the dephasing time - in that case just the 
decay time of Do(t) - is not sensitive to t p . On the other 
hand, a strong coupling regime for which Do(t) decays 
faster than Tlo(t p ,t). In that case, the dephasing time is 
given by the decay time of U (t p , t) and thus possibly in- 
dependent. As the above discussion shows, the cross-over 
between these two regimes may happen for a possibly t p 
dependent cross-over coupling f c (t p ) <C 1. Note that in 
the strong coupling regime, the dephasing time becomes 
independent of the amplitude of the noise whereas in the 
weak coupling regime, it is expected to depend on the 
amplitude and to diverge in the vanishing coupling limit. 



IV. DEPHASING IN THE SYMMETRIC 
MODELS 



In this section we present our results for the situation 
of a symmetric distribution of the spikes, p(— x) = p(x). 
Under this assumption the average random accumulated 
phase vanishes, i.e. f and consequently D t (t) are real. 
For g <C 1, / may be expanded in moments of p(x), 
/ « j 2 « 1 where g measures the typical scale of the 



fluctuations of the spikes. 



A. Dephasing at t p — 

Before discussing the decoherence factor for arbitrary 
preparation time t p , it is useful to investigate the case 
t p = 0. Rewriting the Laplace transform of Do(t) (|27|l as 



(31) 



suggests to introduce a scale 7 g related to the strength 
of the coupling: 



i - L[i>] 



(32) 



s=7s 



Note that 7 9 vanishes with the coupling g. Investigat- 
ing the behavior of Dq for t <^ tq requires evaluating 
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the Laplace transform for stq « 1. In this regime, the 
Laplace transform L[ip]/(l — L[ip]) can be approximated 
by 1/(1 — L[tjj]). Within this approximation, we shall 
now derive explicit expressions for L[-Do] which can be 
Laplace inverted explicitly. This leads to expressions for 
Do(t) valid at t ^> tq for the different classes of /i. 

In the case of finite average waiting time (/i > 1), we 
can expand L[ip] m I — s(r) to find the leading contribu- 
tion to D (t) for (t) < j g t < 1. This gives: 



D (t) 



s(r)+f\ 



,-*/■> 



(33) 



with the dephasing time = 7~ 1 = (t) / / (sec cq. 
(|32(l L Note that this expression is exact only for ip(r) = 
j^j£~ 7 *' ' • At finite non integer /i, taking into account 
the fluctuations of the waiting times requires keeping 
all terms in 1 — L[ip] up to the first non integer power 
{sT b )^. For 1 < [i < 2, we get algebraic subleading cor- 
rections to (j22J|: log(A)(£)) - ~7gt- fc{n)(t/T ) 2 ->*, 
where c(p) = (/t — l)/(2 — /i). These corrections being 
weak for g -C 1, the dephasing time remains equal to 
7~ = (t) I ' g 2 in this regime. 

For jj, < 1, the first term in the expansion of 1 — L[tp] 
is proportional to {t sY . Therefore, plugging in (1 — 
L[i/ji]) / L[i/ji\ ~ r(l — /i)(sro) M and performing the inverse 
Laplace transform of (|3*T|) gives: 



D (t) ~ L- 1 



1 + /(STQ)-" 



= ^[-(7 ff *) M ] (34) 



where E^(z) = Y^=o r(i+^m) denotes the Mittag-LcfHcr 
function^. For values of fi close to one and z <■ 1 this 
function can be approximated by a simple exponential, 
whereas for fi ~ it is similar to an algebraic func- 
tion E^{z) ~ (1 + z) -1 . For large values of the argu- 
ment {\z\ ^> 1) and |arg(— z)\ < (1 — [i/2)ir, we obtain^ 
Efj,(z) « (— zr(l — /a)) -1 . This change of behavior from 
an exponential to an algebraic behavior was interpreted 
in ref. |3^ as a Griffith effect. 

Computing long time behavior {^i g t S> 1) of the de- 
phasing factor can be done by expanding i|31|) as follows: 



Tin V \ ~ 1 1 ~ L W 
L[D \{s) ~ 



(35) 



For 7 s t > 1 (s < 7 S ) we can safely replace « 1 
in the denominator. The inverse Laplace transform can 
then be done easily and reads for j g t ^> 1 : 



1 f°° 
D (t) ~ - J V(r) 



= - 



ro 



f \T +t 



(36) 



For [i < 1 , l|36(l is nothing but the asymptotic behavior of 
(|51Jl for 7 9 t > 1. Note that for g > 1, 7~ x is of the order 
of ro and the decay is algebraic at all times t > tq, given 
by (|SrJ|) . For 1 < [i < 2 and g < 1, only at large times 
("f g t ~ ln(l//)), when the qubit has almost completely 
dephased, the decay crosses over to the power law (I36f) . 



B. Influence of a finite preparation time t v 

We will now discuss in detail the independence of the 
dephasing scenario and of the crossover coupling strength 
g c for the different classes of /i. Simple analytical expres- 
sions for D tp (t) can be derived in the weak (g <C g c ) and 
strong {g ^> g c ) coupling regimes. 



1. Infinite fluctuations of the waiting times: 1 < fj, < 2 

For 1 < /i < 2, the decay of Do(t) is accurately de- 
scribed by and (J35J. On the other hand, U (t p ,t) 
exhibits an explicit dependence on the age of the noise 
(see appendix IP)l : 



U (t p ,t) 



TQ+t 



TO 



T + t + t p 



(37) 



for t p 3> To. We will arbitrarily define the typical time 
scale of the decay of any function as the time at which 
its modulus reaches a fixed value < a < 1 (a = 1/e in 
all the figures of this paper) . The cross-over between a 
weak and strong coupling regime is defined as the value 
of g where the typical decay times of Do(t) and Ylo(t p , t) 
coincide (see section HUB 20 . In the present case, the 
cross-over coupling gets a weak t p dependence: 



1 1/O-1) 



(38) 



Note that g c is a decreasing function of t p since increasing 
tp slows down the decay of Ylo(t p ,t) (remember that the 
average time of the first occurring spike after t p increases 
as ip _A< ). But since the average waiting time is finite, 
g c has a non zero lower bound. Note also that the in- 
dependence of g c is only visible for values of /a close to 
one and disappears as fj, increases to higher values. This 
can be seen on the numerical results depicted on fig. [5] 
g c (t p ) is the crossover coupling where the dephasing time 
start to saturate as a function of g. 

In the case of very strong coupling (g g c (t p )), the 
decay of D tp (t) coincides with Ho(t p ,t) and is thus al- 
gebraic at all times (see fig. |SJ). On the other hand, 
at weak coupling (g <g; g c ), the dephasing factor D tp (t) 
follows Do(t) up to times of the order of the dephasing 
time, i.e. it decays exponentially: D tp (t) = exp(— i/r^), 



where = 7 



Tog . In this weak coupling regime, 



the dephasing is thus described accurately using a second 
cumulant expansion. However, for times large compared 
to the dephasing time, t > t^Ih^t^/tq), higher cumu- 
lants contribute and the decay goes over to a much slower 
power law. If t p < t§ the asymptotic decay of D t (t) con- 
tinues to follow Do(t) behavior given in l|36|l whereas, in 
the opposite case t p > it is given by Ho(t p ,t) l|37|l . 
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FIG. 5: Dephasing factor Dt p (t) in the symmetric model 
obtained by numerical Laplace inversion for fj, = 1.1 and t p = 
10 and various values of g. At small couplings, the decay 
is exponential. But at longer times t S> i> and for strong 
coupling the decay is algebraic. 



FIG. 7: Dephasing time t</, in the symmetric model as a 
function of g for /i < 1 and for various values of t v . At strong 
coupling the dephasing time exhibits a £ p -dependence, which 
disappears as fi increases to higher values. Note the explicit 
tp-dependence of the critical coupling g c . 
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FIG. 6: Dephasing time in the symmetric model as a 
function of g for fi — 1.1, /i = 1.8 and for various values 
of t p . At strong coupling the dephasing time exhibits tp- 
dependence, which disappears as increases to higher values. 
Note that the critical coupling g c has a weak t p -dependence 
for n close to one. 



2. Infinite average waiting time: < fj, < 1 

For < fj, < 1 the influence of a finite preparation time 
t p becomes even more drastic, due to the absence of a 
characteristic time scale in the waiting time distribution. 
In this case ILj(ip, t) only depends on the ratio t/t p (aging 
behavior): 

n»(t P ,i) = ^ 2 F 1 (l,Kl + K-%), (39) 

7T/i \ t J t 

where 2-P1 denotes a hypergeometric functional. Conse- 
quently, the typical decay time of IIo(tp, t) is proportional 
to t p . On the other hand, D (t) is given by (|34f> and its 
decay time thus scales as = jg 1 oc t // 1//ai f» T /g 2 ^ 
in the limit g <C 1. As a consequence, the crossover cou- 



pling strength g c exhibits a strong t p dependence: 

g c (t p ) = A(/i) (f) « 1 , (40) 

where A(/z) is a function of /1 that can be obtained by in- 
version of ffity and As in the case 1 < fx < 2, g c is 
a decreasing function of t p : the range of the strong cou- 
pling regime increases with the age of the noise. However, 
contrarily to the case 1 < /x < 2, g c has no lower bound, 
i.e. it decays to zero as we increase t p . Consequently, 
any qubit surrounded by a noise with < [i < 1 will 
eventually end up in the strong coupling regime. This 
can be seen on results depicted on figure [7| 

At strong coupling (g > g c {t p )), D tp (t) is given by 
Ho(tp, t) and therefore is only a function of t/t p . In this 
regime, the dephasing time is proportional to t p , as shown 
on fig. |H1 The initial decay of coherence is quite fast 
since, for t < t p , U (t p ,t) « 1 - A(^){t/t v ) 1 ^^ with 
A(n) = sin(7r/i)/[(l — /i)7r]. Consequently, for /j, close 
to one Dt (t) decays substantially for times short com- 
pared to the preparation time t p . For t > t p , the decay 
slows down considerably and goes over to a power law 

u (t p ,t) « {tp/ty. 

At weak coupling (g < g c (t p )), the decay of D tp (t) 
is accurately described by (|34|) . For fj, close to 1 and 
7 g £ < 1, the Mittag-Lcfucr function can be approximated 
by a exponential decay D tp (t) « exp[— (7 9 t) M ], whereas 
for 0</i<l and j g t < 1, the decay is rather algebraic, 
D tp (t) ~ [1 + (7 g t)''] _1 . In any case, the typical decay 
time scales as = 7~ 1 oc g~ 2 ^. As for 1 < \x < 2, this 
dephasing time can be recovered using a second cumulant 
expansion of the phase. Obviously, for larger times, t > 
T0, higher cumulants contribute and the decay goes over 
to a power law, D tp (t) oc (ro/t) M , as shown in Fig JHJl. 
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FIG. 8: Dephasing factor Dt p {t) in the symmetric model 
obtained by numerical Laplace inversion for jj, = 0.8, t v = 10 J 
and various values of g. For weak coupling, g < g c (t p ), the 
decay is exponential for t < t$. For times t > at weak 
coupling and for strong coupling the decay is algebraic. 

3. On the marginal case: /i = 1 

Finally, in the marginal case /Lt = 1 corresponding 
to 1/f noise, the above analysis is confirmed qualita- 
tively. We mainly find logarithmic corrections to the 
above results, leading to behaviors intermediate between 
the classes \i < 1 and 1 < \i < 2. The crossover coupling 
g c in this case turns out to be 

/ \ 1/2 

9c(t P )= (^ln(i p /r )J , (41) 

showing a strong dependence on the preparation time t p . 
As in the case fx < 1, it decays to zero for t p — » oo. 

Computing D (t) requires to invert (pTTJl . This can- 
not be done exactly because of the complicated form of 
L[tjj] ~ 1 — stq |ln(sTo)| for stq <C 1. However, Do(t) 
may be estimated considering ln(sro) as a constant and 
replacing s by 1/t in its argument. This leads to: 

Hence, for weak coupling, D t it) decays due to the con- 
tribution of Do(t) on a characteristic time scale ~ 
(ro//)|ln/|. Again, for times large compared to the 
dephasing time, the decay slows down considerably and 
goes over to a power law <|5rj|) . D tp (t) oc t~ x . 

At strong coupling (g > g c (t p )), D tp (t) follows 
Ho(t p ,t), which reads in the present case: 

n fa,t)^ 1 l n fl + ^-) . (43) 

]n(l + tp/To) \ TQ + tJ 

Consequently, the dephasing time exhibits a strong de- 
pendence on the preparation time t p . However, since 
Ho(tp, t) is not a function of t/t p the dephasing time does 
not scale linearly with t p (as for /i < 1) but has a weaker 
t p dependence which depends on the parameter a used 
to define {D tv {r^) = a). 



V. DEPHASING IN THE ASYMMETRIC 
MODELS 

In this section, we consider asymmetric distributions 
p{x) of the pulse amplitudes with a non vanishing average 
h = x. This obviously corresponds to the generic case 
for our intermittent noise, but it also appears naturally 
in the context of decohcrence at optimal points discussed 
in section III CI 

Similar expansions of the exact expression Ij28(l can be 
performed to derive the detailed decay of D tp (t) in the 
case of a finite mean value of the noise, x ^ 0. Never- 
theless, analytic expressions are much harder to obtain 
since in this case / = 1 — e* x is complex. Therefore, in 
this section, we will mainly derive scaling laws of r^, the 
critical coupling g c and present numerical computations 
of D t (t) in various regimes. Note that at strong cou- 
pling, the dephasing is due to single events of the noise, 
i.e. D t (i) ~ Ho(t p ,t). As a consequence, dephasing 
only depends on iP(t), not on the details of p(x). In this 
regime, the results reduce to those previously derived for 
the symmetric case (h = x = 0) in the strong coupling 
regime. Differences with respect to the symmetric case 
only arise in the weak coupling regime on which we shall 
focus in the following. 

Before proceeding along, let us recall the notations Q 
for the moments of p(x) : h = x and g 2 = x 2 . At weak 
coupling, / may be expanded in moments of p(x): f sa 
— ih + g 2 /2 + . . . . As expected for h g 2 , the finite 
mean value does not modify the results of the previous 
section, i.e. the dominant dephasing is the same as in 
the h = case (symmetric noise). However, as h gets of 
the same order as g 2 the decay of D tp (t) and the scaling 
of T0 are modified. Understanding precisely the possibly 
nonlinear effect of a small bias on the dephasing factor 
is related to the fluctuation/dissipation issue in CTRW 
and is out of the scope of the present paper. Therefore 
we shall now focus on the case of huge asymmetry where 
/i> g 2 . In this limiting regime, the pulse dispersion can 
be forgotten and h is the only coupling parameter. Note 
that in this specific variant of the asymmetric model, 
&t p (t) is proportional to the number of events N[t p , t p + 
t] that occur between t p and t p + t and therefore, the 
dephasing factor as a function of h is the characteristic 
function of the probability distribution for N[t p ,t p + t]. 
Contrarily to the case of Poissonian telegraph noised 
where a finite mean value of the noise just adds a global 
phase to D t (t), we will see that strong fluctuations in 
the occurrence times of the spikes induce dephasing even 
for a fixed value of the phase pulses. 



A. The < n < 1 class 

In the < (j, < 1 case, we can use the same method 
as for the derivation of l|34|) to find the scaling law of 
as a function of the coupling strength. The dephasing 
factor in the weak coupling regime decays as \Dt p (t)\ = 
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\Ep(iz)\ where z — — h(t/To) tl /T(l — fi) and Ep denotes 
the Mittag-LefHer function previously used. Using the 
series expansion that defines Ep, it can be shown easily 
that \Ep(iz)\ only depends on z 2 . As a consequence, 
the dephasing time in the weak coupling regime scales as 
oc To/h 1 ^, as shown in Fig.El 

The cross-over between weak and strong coupling 
regime is obtained by comparing this time scale with the 
typical decay time of U (t p7 t), which is proportional to 
t p for < fi < 1. Therefore the critical coupling strength 
scales as h c oc t~^. 




FIG. 9: Dephasing time in the asymmetric model as a 
function of the coupling g for a huge asymmetry /i> g 2 , and 
different values of /i < 1. At weak coupling, h < h c {t p ), the 
dephasing time scales as oc l/h 1 '^, whereas it becomes 
proportional to t p for h > h c (t p ). Note the t p -dependence of 
the critical coupling h c , visualized for [i — 0.4. 

Fig. ^| shows D t (t) for /z = 0.8 and different values 
of the coupling strength from the weak coupling regime 
to slightly above h c . The initial decay, up to the de- 
phasing time, follows the previous exponential decay. 
At larger times the decay crosses over to a power law, 
D t (t) oc l/i M , as in the symmetric case. Between these 
two asymptotic behaviors, beatings may occur, due to 
interferences between different noise histories, as already 
seen in other models^. However, in this reference, oscil- 
lations appear in the strong coupling regime. They arise 
from interferences between noise histories associated with 
different initial conditions. Such a strong initial condition 
dependence is expected at strong coupling for a Poisso- 
nian fluctuator. In our case, oscillations are also present 
at t p = in the diffusion regime and therefore cannot be 
explained by an initial condition dependence. Neverthe- 
less they are still associated with an interference effect 
between noise histories. 

To illustrate this point, let us mention that for D (t), 
a numerical investigation shows that in the situation 
considered here (h ^> g 2 ), these oscillations are only 
present for fi > 1/2. More precisely, the probability dis- 
tribution for the accumulated phase &o(t) can be com- 



puted analytically in the diffusion regime. In this limit 
$o(t) = hN[0, t] can be viewed as a positive real number. 
For < fx < 1, its probability distribution is highly non 
Gaussian and given by 



t j ( t/Tj(h) \ 



(44) 



where r^(h) = r (r(l — ^/h)" 1 ^ and Lp denotes the 
fully asymmetric Levy distribution of index (i. For 
fj, < 1/2, it decays monotonically whereas for \i > 1/2, 
it grows towards a maximum before decreasing rapidcly 
for large values. This maximum can be viewed as a pre- 
cursor of the maximum expected for 1 < \x < 2 close the 
average value ht/(r) of the phase (for 1 < /i < 2, Pqj 
is a truncated Levy distribution whereas for fx > 2 it is 
Gaussian). It is precisely this local maximum that leads 
to oscillations in the modulus of the dephasing factor. 




FIG. 10: Decay of |.D tp (f)| in the asymmetric model for fi = 
0.8, t p = 10 3 ro and different values of the coupling strength 
h. Here /i c (10 3 r ) ~ 10" 2 ' 4 

Note that these oscillations mig ht be related to the 
Griffith effect mentioned in ref. l32l 



B. The intermediate class 1 < fj, < 2 

In the intermediate case of diverging variance, 1 < fi < 
2, it is even harder to find scaling laws analytically. Nu- 
merically, one finds that the decay time of D t (t) in the 
weak coupling regime scales as oc To/h. This leads to 
a critical coupling h c (t p ) equal to the r.h.s of (|38|) . At 
weak coupling, the decay roughly follows an exponential 
exp(— i/r^) up to times of the order of r^. At larger times 
it crosses over to a power law, D tp (t) oc 1/t^. Again, 
at intermediate times and in the weak coupling regime 
Dt (t) shows a transcient regime. Note also that the 
cross-over between the weak and strong coupling regimes 
happens in a much larger (/i-dcpendcnt) range of the cou- 
pling strength than in the symmetric case. 
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VI. DISCUSSION AND CONCLUSION 

To conclude, we have proposed a new model for an 
intermittent classical noise with a 1// power spectrum. 
Within this model, the noise consists in a succession 
of pulses, separated by random waiting times. We 
have shown within this context that the intermittence 
associated with the 1/f characteristics implies a non- 
stationarity of the noise. Non-stationarity effects are 
present in some of the regimes of the dephasing of a two 
level system coupled to this noise as summarized on ta- 
bles [I] and ^ for the symmetric model. In particular, in 
the strong coupling regime, the dephasing factor decays 
algebraically in time, with a characteristic time of decay 
(dephasing time) depending on the age t p of the noise. 
On the other hand, in the weak coupling regime an initial 
exponential decay of this dephasing factor is found. How- 
ever, for low frequency noises and contrarily to the usual 
behavior, this exponential is stretched and the dephasing 
time shows an anomalous dependence on the coupling to 
the noise. After this initial decay, the dephasing factor 
decays algebraically. An important observation is that 
the critical coupling strength separating the weak and 
strong coupling regimes does depend on the age of the 
noise when it is non-stationary. More precisely, it de- 
cays to zero with the age, such that any qubit coupled to 
a non-stationary noise will eventually fall in the strong 
coupling regime. 

One should be careful that the dephasing factor that 
we studied is defined as a configuration average, or en- 
semble average over the noise (see eq. (0). In the 
usual experimental protocols, information on the quan- 
tum statistics of the qubit is collected in a given sample 
in successive runs. This corresponds to a time average 
in a given configuration. These two averages do not co- 
incide in general for non-stationary or aging phenomena. 
Indeed, the non-stationarity of the 1// noise that we con- 
sidered is closely related to the weak ergodicity breaking 
found in similar trap models for glassy materials^. Thus 
one should be careful in interpreting the non-stationary 
dephasing factor that we found in some regimes. We 
nevertheless hope that the questions raised by our re- 
sults might lead to some possible experimental setup to 
better characterize the low frequency noise in mesoscopic 
solids. 



APPENDIX A: USEFUL RESULTS ON THE 
SPRINKLING TIME DISTRIBUTION 

1. Definition 



which states that the spike at t is either the first one, 
or follows a previous spike at time t — t' . In Laplace 
transform, this reads: 



L[S}(s) = 



(s) 



(A2) 



Whenever tp(r) nas a finite mean, S(t) is constant, equal 
to 1/(t) with possible large fluctuations (ji > 1). 



2. Explicit expressions 

In the case of a Poissonian waiting time distribution 
tp(r) = -f e~~ fT , S(t) is a constant equal to l/(r) for all 
times. When tp has algebraic tails at long times, we ex- 
pect this result to be modified since the rate of events is 
expected to decrease with the sampling of the algebraic 
tail of i/j. 



a. Case 1 < fj, < 2 

The small s expansion of L[tp](s) is given by: 

LM(«) ~ 1 - At s - r(l - fx) (stoT 

where A is a numerical constant which depends on the 
small time behavior of tp(r). It is given by A = l/(/x — 
1) for the specific case of (JHJ- Using this form, in the 
vanishing s limit, we have: 



L[ip](s) 



r 1 _m^ {TQsr - 1 _ ATQS+(D{s 2 ) 



Atqs 



A 



(A3) 

Taking the inverse Laplace transform gives the following 
asymptotics for t ^ tq: 



S(t) 



For the specific case of 



A- 



i(t) 



S(t)^(T)-l (l+(^) 



(A4) 



The sprinkling time distribution decreases algebraically 
towards its asymptotic value S'(oo) = 1/(t). This alge- 
braic tail is the signature of the slow decay of tp(r) f° r 
very large times. 



This distribution S(t) is defined as the probability den- 
sity that an event occurs exactly at time (date) t. Hence Case [i = 1 
S(t) satisfies the equation 



In this case, the Laplace transform of tp is given by 
S ® = <«*) + [ dt '^ - 05(f) , (Al) L[m = 1 + STq e ^ Ei(sro) (A5) 
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Exponent fi 


Critical coupling gc(t p ) 


Dephasing time 


y yc 


y <^ yc 


1< n < 2 


9=(*p) 2 = ^[f + 


^> 


M-l" 




T> = T g' 2 


r = r o [l/e+(To/t p )' 1 - 1 ]- 1 /^-i) 


/i = l 






^ = 7^|ln/ s | 


t* = ro^p/ro) 1 - 1 /- 


< /u < 1 


Sc(tp) 2 oc 






r = roff" 2 ^ 


7> oc i p 



TABLE I: Critical coupling and dephasing time in the symmetric model for different waiting time distributions. 



Exponent fi 


.9 > 3c 


5 < 9c 


t < r 


t > ^ 


1< p < 2 


{r + t) ( 




exp(— t/r^,) 


t p < : exp(— t/r^) 




1 In 






1/t 


ln(l + i p /T ) 


«P( / S ln((A)) 


< /i < 1 


t < 7> 


t > 


/i ~ 1 : exp[-(t/T )"] 
^ < 1 . i +(t / T ^)/. 




l-A( M )(t/t p ) 




(i P /i) M 



TABLE II: Summary of the different expressions of the decoherence factor Dt p in the symmetric model. 



where Ei denotes the exponential integral function which 
has the following expansion valid for x > (eq.8.214 in 
ref.ET 



Ei(-x) =C + log(.T) + ^ 



71=1 



(-xl 
klk 



(A6) 



Therefore, keeping only the most singular term in the 
limit stq — ► we get: 

-1 



L[S](s) 



(A7) 



stq log(sro) 

The inverse Laplace transform cannot be found exactly 
but its asymptotic behavior at large times can be esti- 
mated as: 

S(t) ~ — - - . (A8) 

y> T log (t/ t ) 1 ; 

In this case, the sprinkling time distribution decays to 
zero when t — > +oo. This means that events are more 
and more rare. 



c. Case < fi < 1 
In this case, we find 



L[S](s) ~ a ^ 



r(i-M) 



which corresponds to 



S'(t) — t^oo 
also, 

to zero when t — > +oo 



sin(7r/i) 1 /tqX 1- ^ 



(t) 



(A9) 



(A10) 



7T T 

In this case also, the sprinkling time distribution decays 



3. Translated sprinkling time distribution 

Let us denote by St At) the density of events at time 
t p + t. By definition St (t) = S(t p + t). Using the same 
idea as above, we can decompose noise histories in two 
classes: the one that do not have any event between t p 
and t p + t and the ones who do have. This leads to 
an integral equation expressing St (t) in terms of ip and 
ipt defined as the probability distribution for the time 
between t p and the first event occurring after t p . This 
integral equation, analogous to IjAlf) is: 



S tp (t)=^ tp (t) + i>{t-T)S tp {T)dT. (All) 

Jo 

Taking its Laplace transform leads to: 



L[S tp ] 



1-L\p\ 



(A12) 



APPENDIX B: LAPLACE TRANSFORM AND 
MOMENTS OF VM,?/> tp (T) 

1. Laplace Transform and moments of iP(t) 

Let us recall some known results on Laplace Transform 
of algebraic distributions. The distribution (JSJ) can be 
expressed in terms of the reduced variable x = t/tq: 



%jj{x) = fx (1 + x) 



(Bl) 



It is useful to notice that the Laplace transform of (JSJ 
reads exactly 



(B2) 
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Here r(— jx, stq) is the incomplete Gamma function^: 
T(z, a) = f°° e~ t t z ~ 1 dt. The moments of this distri- 
bution exist up to order where [/j] corresponds to the 
largest integer smaller than /i. They read: 

/>oo 

(x n ) =fj, dx ^{l + x}- 1 -^ 
Jo 

= M / dyy n {l-yY- 1 - n 
Jo 

T(n + l)T(u-n) 
= nB(n + 1, /x - n = M 1 > 1 * 

r(/A + i) 

where B(x,y) is the beta function. 

We are interested in the relation between the existence 
of these moments, and the small s behavior of the Laplace 
transform (|B2(I . Expanding the incomplete Gamma func- 
tion r(^, s) in the Laplace Transform (|B2() gives: 



m{ s ) = ^e s r(- M ) 



E 



(-i) fc g fc 



(-l) n - k 



n=0 k=0 



k\(n — k)\(fi + k — n) 



valid for non-intcger \i. The coefficient of s n in this ex- 
pansion reads: 



E 



(— i) n+fe 



=o fc!(n - fc)!(A, + k -n) ^ fc!(n - k)\ J 



E 



(-i) 



(-1) 



n pi 



dt (1 - t)"^" 



= - — f-B(n + l,n-n). 

71! 

Restoring ro, we thus get an explicit expression for the 
Laplace Transform of i/>: 

00 ! — l\ n 

L[1>](s) = M r(-/i)0m ) )'V + £ [ —L( STo r(x n ) (B3) 

The terms n < fi of the last expansion corresponds to 
the moments of ip( x )i an d we recover the expected re- 
sults: the small s expansion of L[ip]{s) consist in a inte- 
ger function that gives the existing moments of tp and a 
second part that is specific of the tails of ip(x) if it decays 
more slowly than an exponential (algebraic tails). 



2. Double Laplace transform of ipt p (Ti) 

In Ref. |H Godreche and Luck used a direct CTRW 
method to derive the double Laplace transform of 
ipt p {i~i). Their results can be recovered straightforwardly 
from the renewal equation (1221) . Wc will use the following 
notation for the double Laplace transform 

L t p ,n[ipt p (Ti)](u,s) = 

dn e - STl / dt p e~ ut * Vt P (n) (B4) 
o Jo 



Let us first focus on the i p -Laplace transform of the inte- 
gral in eq. 1)22(1 . Changing the integration variable from 
t p to t! = t p — T\, it reads 



dtp e~ ut * 



dr S(t p — t)-0(t! + r) 

/>OC 

L[S}{u) / dre-^ipin+r) (B5) 



This provides the following expression for the double 
Laplace transform of ip tp (jl ) : 



L[i> tp \{u,s) = (l + L[S](u)) L tp!T Mtp + ri)Ku,s) 
L t p ,n [ip(t p + n)](u,s) 



(B6) 



The final result can be obtained by considering the 
double Laplace transform of ip(n + t p ) : using 
J^ 00 ip{ T ) e ~ ST dr — L[ip](s) — Jq ip{r)e~ aT dr, we obtain 
(with r = t\ + tp) 

l *p,ti [ip{t p + n)](u,s) 



dr ie -( s -") Tl / dr e-" T V(T) 



«.f, ri r rfTe -Mh e ^ (T) 



s — u 

1 



m](u)-LW\{s)) 



s — u 

Plugging (|B7f) into [|B6|) . we obtain^: 



(B7) 



L Vt p («, s = Lv n VrTTTT B8 

s — u 1 — L[^yJ(u) 

This equation can be used to infer explicit expressions 
for ip t (t) by performing the appropriate inverse Laplace 
transforms (see appendix Q. 



3. Moments of ipt P { T ) 

The expansion of L[ip](s) in powers of s can be used 
to extract the long t p behavior of the moments of ipt (r) . 
For fx > 1, we expect ip tp to have the same algebraic 
decay than tp at infinity. Therefore, only the first [//] 
moments of tpt p arc expected to exist. We will expand the 
double Laplace transform (|B8p of L[ipt ](u, s) into powers 
of s for small values of u. The first coefficients of the 
s expansion corresponds to the Laplace transform with 
respect to t p of the moments of ip tp (t) which we denote 
here by T n (t p ) = (r")^ t . When considering eq. (IB8I) . 
two different contributions from (L[ip] (u) — L[tp] (s))/ (s — 
u) appear. The non integer powers can be expressed as 



E 

fc=0 



E 

k=0 



s ^+k u -k-l 



(B9) 
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while integer powers give a finite sum: 



E 

fe=0 



s k u n - 1 - k . 



(BIO) 



Note that fractional powers are of the form where 
k > 0, which confirms that only the first [fj] moments 
exist. Assuming the following notation: 



• For 1 < /i < 2, the first moment (Ti)(t p ) increases 
sublinearily as obtained from a direct computation. 

• For 2 < /i < 3, the second moment also increases 
sublinearily although it is finite for any t p . Only 
when fi > 3 do we have finite limits for both the 
first and second moments of ipt in the limit t p going 
to infinity. 



(-*)' 



n=0 



(Bll) 



the coefficient of the term in the expansion of 1B8 
in powers of s is given by: 



1 



!£_ u m-i-JV 



E 

m=JV+l 



(-l) m u' 



i-l-JV 



Note that the second term is present only for N < [fx] — l. 
Assuming that fi> 1, 1 — L\ip](u) can safely be replaced 
by (t)u in order to extract the low u behavior of the above 
expression. From this, we infer the Laplace transform of 
T/v(ip): 



L tp [T N ]{v) = {-l) N J^u^- N 



[fj]—m—2 

E 



r(M) (r) 

(-l) m (m + l)! (T w+m ) 
(AT + m)! (r) 



(B12) 



This formula shows that for A^ < [/i], there is a limiting 
value for t v go ing to infinity, given by the w" 1 term in 
the sum (|B12f) . Regular terms that appear in the sum 
(|B12|) contain the short time behavior in t p and, in the 
large t p limit, the non integrer power in the r.h.s gives 
an algebraic decaying contribution. The limiting value of 
7jv(*p) for N < [fi] — 1 is given by: 



APPENDIX C: EXPLICIT EXPRESSIONS FOR 

V>t P (n) 

1. Expression of ipt„(Ti) for 1 < fj, < 2 

The small s expansion of L[i/j](s) reads L[ijj](s) = 1 — 
s(t) +(iT(— J u)(sr ) At (we have to keep all the terms up to 
the first non-integer power). Once again, we set tq = 1 
for simplicity. Plugging it into l|B8() gives: 



L bPt p ](u,s) = (r) -fiT(-fi) 



'1 s^-u^ 1 



s — u J u(t) — iiT{—^i)u^ 
^ (1 + A s"- 1 ) (CI) 



s s(s — u) 



where r has been set to one for simplicity and A = 
M r(-/i)/(r) = -r(l - M )/(r) = r(2 - fi) (we used (r) = 
tq /(/j — 1)). Doing the inverse Laplace transform over s 
yields 



(C2) 



Note that the first constant reads 
At^_ = 1 , 1 / 



(C3) 



r(2 - /i) r V * P , 

Restoring tq, the inverse Laplace transform over u gives: 



Tat(oo) 



(N + 1)(t) 



This result coincides with (|23|l found previously for = 
1 as soon as fi > 1. For A 7 " = [/i], the regular contribution 
to (|B12|) is not there anymore. The [fi\th moment has an 
algebraic sublincar dependance in t p : 



t [m](*p) 



(t) tp 



2. Expression of ipt p (Ti) for < n < 1 



Let us use the small s expansion l|B3(l of L[ip](s) : 
L[ip](s) = 1 + fJ,T(— /i)(sro) At . Using the double Laplace 
(B14) transform l|B8|l and setting tq = 1 for simplicity, we get: 



To summarize, the algebraic tail of ip(r) contaminates 
ipt not only through the divergence of its high moments 
but also through algebraic corrections of the lower ones 
and the sublincar algebraic behavior of the last finite one. 
In particular these results imply: 



^(1-P A 

u — s V V s 



(C5) 



The only assumption in deriving this expression is that 
both t p and T\ are large compared to tq (wc retained 
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only the first terms of the s and u expansions). This 
expression can be exactly Laplace inverted in s : 



(C6) 



We can now perform the inverse Laplace transform over 
u to obtain: 



sin wfx 1 



7T Tl + tp \ T% 



(C7) 



Note that in this case, the short time scale To does not 
appear in if) t (ti) contrarily to the cases fj, > 1. For 
Ti t p , this distribution behaves like 



sin(Tr^) f t P Y n 
V r o/ 



(C8) 



which shows that for very large waiting times, we have 
the same algebraic tail than tp up to a t p dependent renor- 
malization factor. For t\ <C t p , we get another algebraic 
tail : 



■0t p (n) 



sin(7r/i) 1 f t p 



IT t p V Ti 



(C9) 



3. Expression of tp t {r{) for /i = 1 

In the marginal case /n = 1, because of its logarithmic 
variation, S(t p — t) can be replaced by Tolog ^ in the 
integral equation l|22[) . This approximation leads to: 



log(t p /T ) Tl + t p + T Tl + T 



(CIO) 



In the following we always consider situations in which 
tp , ti 3> To , for which the first term in (|C10(I can be 
neglected. Therefore, we get the following asymptotics: 



1 u 



(Cll) 



log(t p /T ) Tl + t p Tl 

Note that in this case, the average (ti)^, 4 is infinite 



APPENDIX D: ANALYTIC RESULTS ON n 

Explicit expressions for Ilo (t p , t) immediately follow 
from its definition and from explicit expressions for ip tp (t) 



in the limit t ^> tq. In the case < fi < 1, the result of 
integration is an hypergeometric function^ : 



IL (t p ,t) 



sin (7r/i) / 1 



TTfi 



j) 2F 1 ( t i,l; 1 



In this case, LIo does not depend anymore on the cutoff 
To and exhibits an aging behavior since it only depends 
on the ratio of time scales t/t p . Using the inversion prop- 
erties of hypergeometric functional (eq. 9.131.1), an al- 
ternative expression can be obtained: 



U (t P ,t) 



sin (vr/i) 



7T/1 



t p + t 



2 F 1 (n,fi;l + fx; 



:) 



tp + t' 
(D2) 

Equation 1D1|) is well suited to the limit t <§C t p whereas 
cq. (|D2|) is better suited to the study of the t p <€. t limit. 

In the case 1 < /j, < 2, the final result still depends on 
To. Using eq. 1|C4() . we get: 



n Q {t p ,t) 



T 



T 



tp + t 



(D3) 



In the limit of vanishing t , Ho(t p ,t) vanishes. Hence 
in this limit, on times scales long compared to f, there 
are always somes events in time intervals of duration t 
contrarily to the case where fi < 1. 

In the /i = 1 case, using cq. <|C11() . we get: 



IL (tt) 



log 1 + f 



log (tp/ro) 



(D4) 



which shares some characteristics with the two previous 
cases. On one hand, as in the 1 < [i < 2 case, it still 
decays to zero at fixed t p and t in the limit To — > 0. On 
the other hand, exactly as in the < yU < 1 case, the 
limits t/tp — * and t/t p — > 1 satisfy: 



lim U (t p ,t) = 0. 

t— >oo 

lim IIo(ip, i) = 1 . 

t v — >oo 



(D5) 
(D6) 
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